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5h ' Abstract 

We examine the behaviour of the polarised structure function data 
for gi(x) in the region of small x (~ 0.01 to 0.1) for the various avail- 
able targets (proton, neutron and deuteron) with the aim of clarifying 
what may be safely deduced with regard to the relation of the cur- 
rently attained small- x region and the asymptotic behaviour as x — > 0. 
We find that fits using a single power-like term are susceptible to the 
isospin combination used. Double-power fits are more stable and also 
provide some evidence for an underlying SU(2) symmetry of the struc- 
ture functions. 
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1 Introduction 



As the data on polarised nucleoli deep-inelastic structure functions improves 
in precision, it becomes ever-more desirable to understand the control that is 
possible over the asymptotic behaviour as x — > for fixed Q 2 . This desire is 
driven by the necessity to extrapolate to the point x = in order to compute 
(or estimate) integrals of the structure functions (as required to test, e.g., the 
Bjorken sum rule). Indeed, while it is generally held that such integrals are 
convergent, a steeply rising behaviour of the type a; -0 ' 9 , for example, would 
leave a significant part of the related integral in the unmeasured region and 
could invalidate any deductions derived therefrom. 

The problem has already been examined |],[2| using varying degrees of 
model-dependent analysis. In this letter we shall attempt to analyse the situ- 
ation starting from the least mo del- dependent approach and then investigate 
the effects of various possible assumptions. In particular, we shall exam- 
ine to what extent simultaneous double-power fits of all the data provide a 
more coherent picture as compared to combinations of single-power fits to the 
single data sets; and what indications there exist for an underlying isospin 
structure. 

We shall avoid any direct reference to Regge theory (save that the struc- 
ture functions should behave as powers of x asymptotically), the effect of 
including model input from Regge-based analyses has been discussed in @. 
We note in passing that in 0] the basic input was a logarithmic behaviour of 
the leading isoscalar contribution. However, given the very short lever-arm 
available in x (the SLAC data between about x = 0.01 and 0.15 were used), 
any logarithmic variation can easily be approximately reproduced by a power 
behaviour and is therefore not to be considered fundamental at this stage of 
the analysis. 

The layout of this letter is as follows: in Section |2] we discuss the sta- 
bility of single-power fits to differing combinations of the data, in Section ^ 
the analysis is extended to include double-power contributions and in Sec- 
tion [| the effect of isospin assumptions is examined. We end with a series of 
conclusions and comments. 



2 Single-Power Fits 



In [|TJ an effectively single-power fit for g™ was performed, with the conclusion 
that the dominant power behaviour was of the form 



-0.015a; 



-0.87 
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If such were the case, the unmeasured region, below x ~ 0.01 say, would 
contribute ~0.06 to the g\ integral. Since even the region below x ~ 0.001 
would still contribute ~ 0.05, this would render any future test of the integral 
on very unsure footing. To save the Bjorken integral from unwanted impli- 
cations, in ]lj the authors were careful to estimate the triplet g\{x) — g™{x) 
combination independently with a smaller, negative, power term ~ —0.45. 

Thus, to provide a starting point, we first perform individual fits to the 
SLAC E143 data at 5 GeV 2 for the three experimental targets: proton, neu- 
tron and deuteron || using a single-power contribution of the form 

9i = ax 5 , (2) 

where a and 5 are then free parameters. The results are summarised in 
Table |], where we also include the hypothetical isovector target p — n. All 
fits return perfectly acceptable values of x 2 although the powers found differ 
widely. A comparison of the fitted curves and the data is presented in Fig. |I| 

Table 1: Results for single-power fits to the SLAC E143 proton, neutron and 
deuteron data, at 5 GeV 2 as described in the text. 



target 


a 


S 


p 
9i 


0.18 ±0.05 


-0.24 ±0.11 


9 n i 


-0.02 ±0.02 


-0.81 ±0.40 


gf 


0.72 ± 1.05 


0.78 ±0.61 




0.10 ±0.05 


-0.51 ±0.18 



As is clear from the results displayed in Table [I], the precise power re- 
turned in a single-power fit is highly dependent on the target combination 
used for the data, varying between —0.81 for the neutron and 0.78 for the 
deuteron. Thus, it is not possible, even in a first approximation to assume 
one universal power to describe this region of x. While this is clearly an 
indication of the necessity to use at least two powers, in view of the results 
of, e.g., [|IJ, where the neutron emerged well described by a single power (as 
too here), it will be instructive to study the variation of the power returned 
as a function of the target considered. To this aim we perform an analysis 
by fitting single powers to combinations of the data according to the formula 

g[ = wi6&-<xxe& (3) 

and vary the parameter 9 from to 180°, thus covering the full range includ- 
ing both gf and the combination g\ — g\. The results are displayed in Fig. |2|. 
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Figure 1: Comparison of the single-power fitted curves and the data. 



The obvious feature is the sharp discontinuity located at 6 ~ 120°, which 
clearly indicates the failure of a single-power fit in that region of combina- 
tions, very close to the deuteron. What should also not be underestimated 
is the tail of this discontinuity, which extends into the region of pure g™. 

As commented near the end of the paper, the few (2 or 3) smallest-x 
data points (below x ~ 0.03) of the E155 Q (and SMC) experiments (at 
5 GeV 2 ) hint that more small- x data in the future may push the deuteron 
fitted curves towards negative powers, which on the 9 plot of a single-power 
fit would mean a discontinuity and thus inconsistency. 

A consistent set of fits, related to the procedure of 0, is provided by 
taking the functions fitted to g\ and g\ — g™ as "primary" and obtaining 
g\ and gf as linear combinations, so that there are only two independent 
powers (all functions have acceptable % 2 ). As the starting point may be any 
two of the four functions (with single powers), there is to be an element of 
arbitrariness here. Such arbitrariness is implicitly exploited in |TJ where the 
single-power fit to g± is considered more fundamental or "primary", thus 
leading to a very divergent g\. However, such behaviour disappears if g\ is 
taken as a linear combination, as also intuitively suggested in ||. Is there 
anything then that favours two of the functions and their associated single- 
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mixing angle, 6 

Figure 2: The power, 5, obtained in a single-power fit as a function of 8, the 
mixing parameter for g\ and g™ described in the text. 

powers as more fundamental than the other two? 

Here, we are led once again to consider two separate power-term contribu- 
tions, at least in two of the four combination gi functions considered. Thus, 
it is reasonable to begin with a two-power global fit to the structure func- 
tions to avoid all arbitrariness. This ought to be done simultaneously for all 
the available data. We actually used SLAC and HERMES at approximately 
3GeV 2 . 

Before moving on to double-power fits, it is worth investigating the pos- 
sibility of describing both the proton and neutron data with the same single 
power {i.e., explicitly leaving aside the deuteron data for the moment); the 
results are 

a p = 0.10 ±0.02 

a n = -0.05 ±0.01 (4) 
5 = -0.48 ±0.08. 

However, the \ 2 here is significantly poorer than in any other fit so far. On 
examining these results together with the data, it emerges that the poor fit 
could be readily remedied by the inclusion of a plateau increasing the overall 
proton polarisation and decreasing that of the neutron. The effect would then 
be to reduce the neutron power while increasing that of the proton, leaving 
a slightly positive deuteron (given by the plateau contribution itself). 
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3 Double-Power Fits 



Not only is it obviously necessary to allow for more than one power in the 
fits, but it should also be clear that any single power returned above could 
change dramatically in a two-power fit ||. In particular, two facts lead to 
the possibility that any given behaviour might be mimicked by combinations 
of different powers. 

First of all is the limited lever-arm available in terms of the x-range 
of the data one can sensibly use; little more than one order of magnitude 
variation in x is the bare minimum for extracting a power behaviour, given 
the present experimental errors. Second, and often neglected, is the fact 
that, in contrast to the unpolarised case, where most of the standard wisdom 
has been gathered, the polarised structure functions are not constrained to 
be positive definite. This means that, e.g., an overall rising behaviour (in 
magnitude) may be due either to a single rising contribution or to a rising 
cancellation between two relatively changing contributions of opposite sign, 
neither of which need necessarily be similarly rising. 

In H an explicit example was invented ad hoc, in which a behaviour of 
the type (as proposed in 

n, ^ -°- 02 /k\ 

9i( x ) ~ ~^nr ( 5 ) 

was shown to be well reproduced, within experimental errors and over the 
finite range of x considered, with the form 

which, one should note, has very different asymptotic behaviour. 

The fact that single-power fits do, in fact, work rather well implies that 
an attempt to fit any single target data set (excepting perhaps that of the 
deuteron) with two power terms will encounter serious difficulties. Indeed, 
the only way to perform successful two-power fits is to combine data sets of 
different targets and demand that the two powers used be the same. The 
results of such a fit can be summarised as follows: for the form 

gx = ax 5 + fix*, (7) 

and the SLAC M E143, E154 together with the HERMES M data sets, 
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we obtain 



Op 

5 

Pn 

7 



0.01 ±0.02 
-0.02 ±0.03 
-0.77 ±0.33 
0.26 ±0.20 
0.00 ±0.10 
0.13 ±0.50 



and again the y 2 is perfectly acceptable. Note that the large errors are mainly 
due to strong correlations between coefficients and powers and therefore ul- 
timately between the coefficients themselves. Again, a comparison of the fits 
and data is shown in Fig. 0. 
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Figure 3: Comparison of the double-power fitted curves (solid lines) and 
isospin type fits (dashed lines) with the data. 



One sees the attraction of the neutron data to a single-power fit, i.e., that 
a second power is not at all required by this particular target. Moreover, the 
fact that the power chosen is very close to the original power but far from 
that of the proton is indicative of the possibility to mimic powers (in this 
case that of the proton) by combining different contributions. Note also that, 
as has been found by various authors, the steeper behaviour (5 here) , having 
opposite sign coefficients (a) in the proton and neutron, is to be ascribed to 
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an isovector contribution — within errors the two coefficients have the same 
magnitude. 

As a final test we fix the steep power to be less divergent, setting 5 = —0.5, 
and repeat the fits — the idea being to present a picture less clouded by the 
correlations between parameters, which often artificially inflate errors; the 
results are 

a p = 0.03 ±0.05 
a n = -0.07 ±0.03 

5 = -0.50 fixed . . 

P P = 0.17 ±0.09 W 
(3 n = 0.11 ±0.07 
7 = 0.02 ±0.63 

as always the \ 2 is perfectly acceptable and indeed marginally improves due 
to the reduced number of parameters. 

Hence it is clear that such steep powers as previously found are not an 
absolute requirement of the data, as suggested by the large errors found ear- 
lier. Moreover, on choosing a less divergent behaviour the various coefficients 
come more into line with isospin symmetry requirements; i.e., one term is 
approximately isovector and the other isoscalar. 



4 Isospin Symmetric Fits 

In the light of the above fits, it is natural to enquire as to the effect of 
requiring that the two terms used fall precisely into the two categories of 
isovector and isoscalar. Thus, as a final test we have simply repeated the 
above fit fixing a n = —a p and j3 n = (3 p . The results are (with the usual good 

x 2 ) 

a p = 0.78 ±0.02 

6 = -0.45 ±0.07 , , 
(3 P = 0.20 ±0.14 1 ' 

7 = 0.36 ±0.28, 

where 5 and 7 now refer to the isovector and isoscalar powers respectively. 

Relaxing one or other of the isospin constraints leads to similar fits with 
similar \ 2 '■ Note that in this case the steep power, which is, of course, 
now constrained to be that of the combination g\ — is much less steep 
than in many other fits. This would mean an increasing difference in total 
polarisation of u and d quarks, with decreasing x in the given interval. A 
comparison of the fits and data can be found in Fig. [| 

Thus, we may say that in the small-x region under examination, the 
double-power fits of the data prior to the E155 experiments are consistent 
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with the assumption that the positive isotriplet contribution is causing the 
mildly divergent behaviour of gi(x) and g™ (x). However, the latest E155 data 
on the deutron M, which are still preliminary, provide a hint of a negative 
singlet term in gf(x) with a more strongly divergent behaviour, setting in 
below x = 0.03, in accordance with the expected || asymptotic dominance 
of the singlet. This, together with the observed attraction of the neutron 
data to a single steep power, calls for a fresh look at the picture when the 
final E155 data become available. 

5 Conclusions 

The first lesson that we wish to underline is the possibility that any result 
arising from an effectively single-power fit may be biased precisely by the 
choice of a single term to describe the data. Thus, although the neutron data 
appears to select a large single power, two lesser powers are also perfectly 
acceptable (and less dramatic). Indeed, by comparing the magnitudes of the 
isoscalar and isovector contributions, one sees that over much of the x range 
of interest neither contribution is negligible. It is also clear that the presence 
of an isoscalar plateau leads to a steepening of one power and a flattening of 
the other. 

Perhaps, the second lesson is that, while the isospin properties of the 
distributions are still far from well defined in the data, a simple description 
based on such a symmetry does, in fact, work rather well. 

One caveat that should be borne in mind is that we have not considered 
how evolution might alter the picture: this is both beyond the scope of the 
present letter and, to any reliable degree, beyond the level of the present 
data given the poor determination of the gluon distribution. 
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